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ABSTRACT 

Confexf. Energetic feedback from Supernovae and stellar winds can drive galactic winds. Dwarf galaxies, due to their shallower 
potential wells, are assumed to be more vulnerable to this phenomenon. Metal loss through galactic winds is also commonly invoked 
to explain the low metal content of dwarf galaxies. 

Aims. Our main aim in this paper is to show that galactic mass cannot be the only parameter determining the fraction of metals lost 
by a galaxy. In particular, the distribution of gas must play an equally important role. 

Methods. We perform 2-D chemo-dynamical simulations of galaxies characterized by different gas distributions, masses and gas 
fractions. 

Results. The gas distribution can change the fraction of lost metals through galactic winds by up to one order of magnitude. In 
particular, disk-like galaxies tend to loose metals more easily than roundish ones. Consequently, also the final metallicities attained 
by models with the same mass but with different gas distributions can vary by up to one dex. Confirming previous studies, we also 
show that the fate of gas and freshly produced metals strongly depends on the mass of the galaxy. Smaller galaxies (with shallower 
potential wells) more easily develop large-scale outflows, therefore the fraction of lost metals tends to be higher. 

Key words. Galaxies: abundances - Galaxies: dwarf - Galaxies: evolution - Galaxies: ISM - Galaxies: jets 



1. Introduction 

Theories of cold dark matter-dominated hierarchical growth of 
structures predict that dwarf galaxy- (DG-)sized objects are the 
building blocks for the formation of large galaxies. In spite of 
their relevance, the most important physical phenomena regulat- 
ing the birth and evolution of DGs are still obscure to date. For 
sure, star formation (SF) plays a key role in shaping DGs and 
determining their fates. Since the binding energy of the interstel- 
lar medium (ISM) in early gas-rich DGs is very small (smaller 
than the explosion energy of just a few Supernovae (SNe)) many 
authors speculated that a high SF rate in a DG would create a 
galactic wind and thus produce a transition from a dwarf irreg- 
ular (dlrr) to a dwarf spheroidal (dSph) or dwarf elliptical (dE) 
■ (see e.g. Larson [19741 Vader [T986l Dekel & Silk fTWSl From a 
chemical point of view, the occurrence of a galactic wind right 
after the formation of the first stars would imply a very lim- 
ited interval of time during which the metals, restored by dy- 
ing stars, can pollute the ISM and enrich the following stellar 
populations. If galactic winds preferentially occur in DGs, these 
objects (with low masses) will experience a limited chemical en- 
richment. Therefore, a correlation between stellar mass M (or 
luminosity L) and metallicity Z of DGs is expected. Indeed, this 
M-Z correlation among DGs exists (see e.g. Skillman et al. 1989 1 
Lee et al. 2006; Kirby et al. 120081 Zhao et al. 120101 Andrews & 
Martini l2012l ). A M-Z relation extends also to high redshifts (Erb 
et al.120061 Maiolino et al. l2008l Laskar et al. BOTTl Wuyts et al. 
120121 although, in this case, the galactic masses of the galax- 
ies for which metallicity determinations are available are usu- 
ally quite high (but see Mannucci et al. 1201 II) . The M-Z relation 
among DGs corroborates the idea that SN-driven galactic winds 
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play a dominant role in the evolution of these objects. More re- 
cently, a general relation between stellar mass, gas-phase metal- 
licity and star formation rate has been found (Mannucci et al. 
l2TJT0l l. 

However, detailed hydrodynamical simulations of DGs 
showed that the galactic winds, although often able to expel a 
large fraction of freshly produced metals, are unable to eject an 
equally large fraction of pristine (i.e. not processed) gas. This 
is mostly due to the fact that, if the initial DG gas distribu- 
tion is flattened (as observed in dlrrs), then there is a direc- 
tion with steeper pressure gradient so that the galactic wind will 
preferentially expand along that direction, the transport of gas 
along the other directions being very limited (see e.g. D'Ercole 
& Brighenti [19991 MacLow & Ferrara IT9991 (hereafter MF99); 
Recchi et al. 2001). This effect can be appreciated by inspect- 
ing Fig. Q] where the gas and oxygen distribution of a model 
galaxy experiencing a galactic wind (taken from the calculations 
of Recchi et al. 120061 ) is shown. Most of the disk gas (at least 
above R= 1 kpc) has not been affected by the galactic wind and 
it will fall back towards the center of the galaxy once the energy 
source of the starburst will be exhausted. On the other hand, the 
oxygen can be easily channelled along the funnel created by the 
galactic wind. Also very energetic (E ~ 10 53 erg) hypernovae 
exploding in dwarf protogalaxies, although more disruptive than 
normal SNe, are able to expel only ~ 10 to 20 % of the baryonic 
mass originally present in the galaxy (Vasiliev et al. 2008 ). 

This result is widely accepted by the astrophysical commu- 
nity, although some authors (e.g. Tenorio-Tagle 1996; Silich & 
Tenorio-Tagle 1998 ), believe that low-density gaseous galactic 
halos can strongly affect the circulation of the metal-rich mat- 
ter processed by the central starburst. Kiloparsec-scale bipolar 
outflows are created but, later on, loss of pressure support and 
interaction with the diffuse halo slow down the expansion of the 
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Fig. 1. Gas (filled contours) and oxygen (black contour lines) 
distribution of a model galaxy (taken from Recchi et al. 2006) 
experiencing a galactic outflow. Brighter colors indicate larger 
gas densities. 



superbubble and the freshly produced metals can eventually fall 
back towards the center of the galaxy. Thus, according to this 
scenario, energetic events associated to a starburst can still cre- 
ate a large-scale outflow but this outflow does not necessarily 
become a galactic wind (i.e. its velocity does not exceed the es- 
cape velocity of the galaxy). 

It is also worth noticing that, while the chemical yields of 
exploding massive stars are very large (thus the metallicity of 
interiors of SN-driven bubbles is up to 40 times the solar value), 
abundances of the hot gas in galaxies experiencing large-scale 
outflows (determined from X-ray spectra) cover with high proba- 
bility a range between only solar and twice solar (see e.g. Martin 
et al. [2002for NGC1569 or Ott et al. [2005]for a larger sample 
of galaxies). From this fact, it can be deduced that the hot SN 
ejecta becomes mass loaded (by a factor of up to 10) during 
its expansion, probably due to evaporative and turbulent mix- 
ing with colder interstellar clouds. By this, the hot gas loses not 
only momentum but also energy by enhanced cooling in addi- 
tion to adiabatic expansion so that the outbreak of the hot gas to 
a galactic wind could be hampered in many cases. This mass- 
loading effect is probably causing the discrepancy between the 
observed extent and temperature of the extremely cool superbub- 
ble in NGC 1705 and its analytic parameter correlation (Hensler 
et al. [19981 

Because of the gas density outside the SF region, the de- 
velopment of galactic outflows or galactic winds considerably 
depends on the degree of flattening (or on the rotation) of the 
parent galaxy, an aspect that has not been fully explored in the 
past (but see Silich & Tenorio-Tagle 2001, Michielsen et al. 
120071 Schroyen et al.l20TT1 See also S trickla nd & Steven s 120001 
Ferrara & Tolstoy |2000] Vasiliev et al. 120081 Recchi et al. 120091 . 
Although flat galaxies are supposed to loose a large fraction 
of metals as a consequence of bipolar galactic winds, roundish 
galaxies are not characterized by a direction along which the 
pressure gradient is significantly steeper, therefore they are more 
likely to retain much of the metals (Marcolini et al. 2006 , Recchi 
et al. 120071 1 . In the end, this difference in the initial gas distribu- 
tion could lead to a spread in the observed final metallicity of 
galaxies with the same mass. Such a spread is indeed observed 



in the M-Z relation of DGs (Lee et al. 120061 Kirby et al.[2008; 
Zhao et al. 120101 . or in the same correlation at z — 3.5 (Maiolino 
et al. I2008I I. In this paper we explore in detail the effect of gas 
geometry on the development of galactic winds and on the fate 
of freshly produced metals. This work must be thus seen as a re- 
fined study along the lines of the MF99 work, in which variations 
in the initial distribution of gas in galaxies have not been consid- 
ered. The second key parameters we consider in this study is the 
initial mass of the galaxy. As already mentioned, small galaxies 
(with shallow potential wells), are expected to develop galactic 
winds more easily than larger galaxies, therefore the mass is cer- 
tainly a key parameter in determining the fate of pristine gas and 
freshly produced heavy elements. 

The dynamical evolution of a galaxy is so complex that the 
fate of heavy elements cannot depend only on mass and geom- 
etry of the parent galaxy. An obvious parameter that strongly 
affects the development of galactic winds is the luminosity 
(MF99). In turn, the luminosity depends on the star formation 
rate (SFR) and on the star formation history (SFH) and it is 
known, at least in the Local Group, that very large galaxy-to- 
galaxy SFH variations exist (see e.g. Monelli et al. 2010a, 2010b; 
Hidalgo et al. 1201 11 1. Moreover, other details of the ISM struc- 
ture (for instance the presence of clouds or porosity) and of the 
feedback prescriptions (e.g. how the SN explosion energy is re- 
distributed, where, on which timescales and at which efficiency 
etc.) can all play a role in determining the amount of metals car- 
ried out of a galaxy by galactic winds and eventually in the final 
global metallicity of the galaxy. Although the consequences of 
an inhomogeneous ISM on the development of galactic winds is 
analyzed in detail in Recchi & Hensler (2007 hereafter RH07), 
we have undertaken a comprehensive parameter study through 
detailed chemo-dynamical simulations, the results of which will 
be presented elsewhere. Here we will focus mainly on 18 ba- 
sic models, characterized by different geometries, total baryonic 
masses and initial gas fractions. The paper is organized as fol- 
lows: in Sect. [2] a short summary on the literature results con- 
cerning the fate of freshly produced metals will be given. In 
Sect. [3] our numerical scheme will be briefly reviewed and the 
set-up of the 18 basic models models described. In Sect. [4] the 
main results of these calculations will be presented, underlying 
also the effect of some other key parameters not fully explored 
in this work. In Sect. [5] the results will be discussed and some 
conclusions will be drawn. 

2. A brief overview of literature results 

As already mentioned in the Introduction, several papers in the 
literature attempted to study the effect of galactic winds on the 
circulation and redistribution of metals in DGs. The main results 
of the often-cited work of MF99 are that, even in the presence 
of a strong galactic wind driven by SNell, the ejection efficiency 
of unprocessed gas is always close to zero (with the exception 
of galaxies with initial baryonic mass < 10 6 M ). On the other 
hand, the ejection efficiency of freshly produced heavy elements 
is always close to one (with the exception of galaxies with initial 
baryonic mass > 10 9 M ). D'Ercole & Brighenti (1999) found 
results similar to the ones of MF99, in the sense that only mi- 
nor fractions of ISM (but in some cases large fractions of met- 
als) are expelled through galactic winds. They however calcu- 
lated the evolution of the galaxy for much longer times than 
MF99, discovering that the center of the galaxy can be replen- 
ished with cold gas in a timescale of ~ 100 Myr (see also Recchi 
& Hensler 2006). As already mentioned, Silich & Tenorio-Tagle 
(1 19981 1 found instead that in most of the models galactic winds 
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do not develop (mainly due to the presence of a hot gaseous halo 
surrounding the galaxy). Fragile et al. (2004 1 studied the effect of 
SN explosions not localized in the center of the galaxy (as MF99 
did) but distributed over radii of up to 80% of the disk radius, 
discovering that, in this case, radiative losses are more effec- 
tive and the development of a large-scale outflow is hampered. 
Also very interesting is the work of Scannapieco & Briiggen 
( 120101 1. who attempted to model turbulent velocities and turbu- 
lent length scales in DGs and to inject SN energy into supersonic 
turbulence. The wind efficiencies they found are still quite low. 
Rodriguez-Gonzalez et al. (1201 11 1 addressed the same problem 
(mass and metal ejection efficiencies in DGs) exploring a much 
more extended set of parameters (in particular they covered a 
wide range of starburst and galactic gas masses). Their results 
are overall in agreement with the previously quoted works (with 
the exception of Silich & Tenorio-Tagle 1998). More recently, 
Schroyen et al. (201 1) explored in detail the effect of rotation on 
the evolution of DGs. Their focus was to understand the origin 
of the dichotomy in radial stellar metallicity profiles of DGs and 
their "centrifugal barrier mechanism" help explaining these ob- 
servations. Usually their models do not develop galactic winds 
(or develop only weak outflows), thus the determination of ejec- 
tion efficiencies (of pristine gas and freshly produced metals) is 
not addressed in their paper. The same can be said about the de- 
tailed determination of the energy required to expel newly pro- 
cessed matter by Silich & Tenorio-Tagle (2001); their analytic 
model is in good agreement with the numerical results of MF99 
and shows clearly the effect of geometry on the development of 
galactic winds and on the fate of metals if the intergalactic pres- 
sure is low. However, their inferred large intergalactic pressure 
hampers the development of galactic winds. 

As this short summary indicates, the past literature mostly 
focused on the effect of mass and luminosity on the gas and 
metal ejection efficiencies from galaxies. The chemical evolu- 
tion is usually not very detailed in these models. Recchi et al. 
( T200T1 120021 120041 120061 ) provided detailed information on the 
ejection efficiencies of single chemical species produced by SNe 
(of both Type II and Type la) and by winds from intermediate- 
mass stars. RH07 added important information on the effect of 
clouds. The most significant results of these works can be sum- 
marized as follows: 

- Galactic-scale outflows carry out of the galaxy mostly the 
chemical elements produced by dying stars during the most 
recent episodes of SF, with large escape fraction of metals 
with delayed production (like Fe and N). 

- Models with very short burst(s) of SF can cool and mix 
the newly formed metals in a very short timescale, whereas, 
when the SFH is continuous and long-lasting, most of the 
metals are either directly ejected outside the galaxy through 
galactic winds or are confined in a hot medium, therefore 
cannot contribute to the chemical enrichment of the warm 
ionized medium observed by emission lines from the HII 
gas. 

- Models with complex and long-lasting SF episodes taken 
from the literature reproduce the chemical composition and 
the abundance ratios, observed in galaxies, much better than 
models with bursting SF. 

- Clouds (if present) are able to increase the main density of 
the cavity created by the ongoing SF activity, provoking a re- 
duction of the total thermal energy by ~ 20 - 40% compared 
to models without clouds. 

- The interaction clouds-supershell leads to strong structuring 
and piercing of the shell, allowing the venting out of metals 



in spite of the reduced thermal energy. The development of 
large-scale outflows is therefore generally delayed but the 
ejection efficiency of metals remains unchanged. 
- Since the clouds are assumed to have a low metallicity, then- 
mixing with the ISM tends to reduces the abundance of 
heavy elements of the galaxies. The models with infalling 
clouds have thus metallicities up to ~ 0.4 dex lower than 
corresponding models without clouds. The abundance ratios 
remain unaltered. 

A study of the early evolution of tidal DGs (which are as- 
sumed to be dark-matter- (DM-)free) and their possible correla- 
tion with the satellite DGs surrounding the Milky Way has been 
also performed (Recchi et al. 2007). These models, in spite of be- 
ing DM-free, can sustain the energy released by dying stars with- 
out experiencing a complete blow-away, for several hundreds of 
Myr, provided that they keep their initial spherical symmetry. A 
typical model develops through a network of cavities and fila- 
ments due to the patchy distribution of the SF sites. A super- 
shell is formed, but it grows slowly in size and does not quench 
the SF process. This is because the initial ISM distribution in 
the galaxy is spherical and there is therefore no preferential di- 
rection through which the galactic wind can flow. Under these 
conditions, either the feedback is able to remove all the gas at 
once (blow-away) or all the gas (or most of it) is retained inside 
the galaxy. Our study, and other similar ones (e.g. Marcolini et 
al. [2006), show that it is not easy to get rid of all the gas in an 
initially spherical galaxy. 

3. Description of the model 

This work aims at studying in detail the dynamical and chemical 
evolution of DGs and the effect of key parameters (in particu- 
lar baryonic mass, initial gas mass fraction and degree of flat- 
tening) on the development of galactic winds and on the fate 
of freshly produced metals. A follow-up study will extend the 
parameter space and will concentrate on the details of the SFH 
and on the effect of boundary conditions. The ultimate goal is 
to understand the main mechanisms determining the metallicity 
and the abundance ratios in DGs and, at the same time, infer the 
degree of metal pollution of the inter-galactic and intra-cluster 
medium due to these objects. 

3. 1 . The chemo-dynamical code 

As chemo-dynamical code we use basically the one of Recchi 
et al. (12007ft . We recall here the basic features of this code. It is 
a 2-D code in cylindrical coordinates based on a second-order, 
MUSCL-type upwind scheme (the 1-D version of this scheme 
is described in Bedogni & D'Ercole 1986). The chemical en- 
richment of the galaxy is followed in detail; the production (by 
SNe of Type II, Type la and intermediate-mass stars) of 8 chem- 
ical elements (H, He, C, N, O, Mg, Si and Fe) is considered 
and the advection of these elements is followed by means of 
passive scalar fields. For instantaneous bursts, this scheme is 
described in detail in Recchi et al. (120011 1 and its extension to 
more complex SFHs is explained in Recchi et al. ( 120041) . Since 
the code keeps correctly track of the evolution of the metallic- 
ity in each computational cell, the detailed metallicity-dependent 
cooling function of Bohringer & Hensler ( 1989) can be imple- 
mented. The heat transport equation is also solved by means 
of the Crank-Nicolson method (see D'Ercole & Brighenti 1999 
for details) using the classical Spitzer-Harm thermal conductiv- 
ity (Spitzer & Harm 1 19531 Spitzer [T956l . A saturated heat flux 
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(Cowie & McKee 119771) is adopted if the mean free path of 
electrons is larger than the temperature scalelength. Some im- 
provements of this code are described in Recchi et al. (2007i. 
In particular, metallicity-dependent stellar winds from massive 
and intermediate-mass stars have been considered and, primary 
and secondary nucleosynthetic production by stars of any age 
and any metallicity has been carefully calculated. Recchi et al. 
( 120071 1 describe also an implementation of SF recipes. The self- 
consistent SF module is not adopted in this work where, in anal- 
ogy with many previous studies (MF99 for instance) the SFH 
(or the luminosity) is an input of the model and not the result 
of the galactic evolution. In a work in preparation, models with 
a self-consistent SFH will be shown and discussed. Moreover, 
although in Recchi et al. ( 120071 1 a self-gravity solver (based on 
Rieschick & Hensler 2003) has been implemented, we decided 
not to adopt it in this study neither. The main reason for this 
choice is that we neglect self-gravity in building the initial equi- 
librium configuration (see below in the next Subsection). A par- 
allel line of research (Vorobyov et al. 2012, Vorobyov et al., in 
preparation) analyzes in detail DG models with a consistent im- 
plementation of the gas self-gravity. The inclusion of self-gravity 
in the code would imply that, at the beginning of the simula- 
tion, the gas is initially out of equilibrium even without sinks 
or sources of energy. This would establish an inward gas flow. 
Since our main aim in this paper is to study the effects of ge- 
ometry and other parameters on the flow rate of gas and metals, 
we do not want gas infall to affect our results and we want flows 
of gas (and metals) to be solely affected by feedback processes. 
Moreover it is important to notice that most of similar papers 
already described in Sect. |2](e.g. MF99, D'Ercole & Brighenti 
1999 Strickland & Stevens 2000) also neglect self-gravity. 

3.2. The set-up 

We start with a reliable set-up of a DG, in which the gas is ini- 
tially in isothermal equilibrium with a spherical DM halo and 
with the centrifugal force. Analogously to what is done in many 
similar studies, we neglect self-gravity in building the initial 
equilibrium configuration. We outline here that this is in gen- 
eral not correct even if most of the mass of the model galaxy 
is in the form of DM and the correct derivation of an equilib- 
rium configuration is the one outlined by Vorobyov et al. (I2012I I. 
However, the equilibrium configurations produced by Vorobyov 
et al. (120121 ) are quite complex (and also computationally de- 
manding) and it is difficult to obtain two analogous equilibrium 
configurations, differing only for the degree of flattening of the 
initial gas distribution. 

We consider three possible values for the initial baryonic 
mass of the galaxy (10 7 , 10 s and 10 9 M ) with a factor of ~ 10 
more massive DM halos. The exact factor is deduced by the cor- 
relation between the dark matter-to-baryon ratio <p and the gas 
content as adopted by MF99, namely <p - 34.7 M~°- 29 , where 

M g j is the gaseous mass of the galaxy in units of 10 7 M Q ). 
This correlation is adapted from the work of Persic et al. ([1996), 
and is based on fitting procedures of observed rotation curves 
in local galaxies. One should warn the reader that it is not clear 
whether present-day dark matter-to-baryon ratios also correctly 
reflect the initial conditions of the galaxy, but we nevertheless 
adopt this correlation in order to facilitate the comparison with 
the results of MF99. The virial radius of the DM halo is assumed 
to be: 

r vir = 0.75 Ml /3 h' 1 i^-^j kpc, (1) 



(Madau et al. l200Tl Mori et al. 120021 . In this formula, M 8 is the 
halo mass in units of 10 8 M and z g f is the redshift of galaxy for- 
mation. Our reference value for z g f is 8. It is important also to 
stress that our initial distributions of gas are not artificially trun- 
cated at some cut-off radii as the models of MF99. Therefore, 
without sources or sinks of energy, our models preserve their ini- 
tial configurations (at variance with the models of MF99 which 
tend to expand). However, since our distribution of gas extends 
until the edges of the computational box, we need to establish 
a radius within which we calculate the mass of the galaxy. We 
take this to be half of the virial radius. Therefore, galaxies with 
the same nominal baryonic mass are normalized in such a way 
that the total mass within 0.5 • r,,,> is the same, irrespective of the 
degree of flattening. This implies that flat galaxies have larger 
central densities than roundish ones (see below and Table [2}. 

A small stellar disk is also initially present in the galaxy. Its 
density is such that a Miyamoto-Nagai potential 



y]R 2 + (a + V?T/7) 

(where Md is the mass of the stellar disk and R, z are the cylin- 
drical radial and vertical coordinates, respectively), is repro- 
duced. The ratio b/a between the scale lengths identifying the 
Miyamoto-Nagai potential is taken as one of the key parame- 
ters that we vary in our models. A small b/a corresponds to a 
flat model (for b/a — > the potential tends to the razor-thin 
Kuzmin model), whereas if b/a is very large, the galaxy tends to 
be rounder (for b/a — > oo the potential tends to the Plummer's 
spherical potential). 

It is expected that an initially flat distribution of gas results 
in an easier development of bipolar outflows, since the pressure 
gradient along the polar direction is much steeper than that along 
the galactic disk. As already discussed, along the polar direction 
a significant fraction of the matter processed in the central SF re- 
gion can be thus channelled and eventually lost from the galaxy. 
On the other hand, if the galaxy is initially spherical (or almost 
spherical), no preferential propagation direction for the super- 
bubble exists and the freshly produced metals are more likely to 
remain confined inside the galaxy. One of the aims of this work 
is to confirm this empirical assumption and quantitatively deter- 
mine the fraction of gas and metal-rich stellar ejecta lost from a 
DG, as a function of its degree of flattening. We consider three 
representative values for the ratio b/a: 0.2 (flat models, desig- 
nated with the letter "F"); 1 (medium models or "M") and 5 
(roundish models or "R"). However, the mass of this pre-existing 
stellar disk cannot be established a priori, therefore as another 
parameter we vary the ratio between the mass of the Miyamoto- 
Nagai stellar disk and the total baryonic mass initially present in 
the galaxy. In particular, we consider two basics sets of models; 
for the first set (designated with "H") the initial gas fraction is 
high (90% of the total baryonic mass of the galaxy) and, conse- 
quently, the pre-existing stellar disk represents a small fraction 
of the mass budget in the galaxy. The second set of models (des- 
ignated with "L") is characterized by a much smaller initial gas 
fraction (60% of the baryonic mass). Therefore, for instance, the 
model H7M represents a galaxy with high initial gas fraction, 
10 7 M of initial baryonic mass and a medium degree of flat- 
tening (b/a =1). We consider thus a total of 18 basic models, 
but we will discuss the dependence of our results on some key 
parameters in Sect. 14.21 

Other key parameters to be varied are the duration and in- 
tensity of the SF episode (as already mentioned, the analysis of 
models with self-consistent SFHs will be deferred to a future 
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work). We assume here for simplicity that the SF is constant 
(with some pre-defined intensity if/) for a period of time At. The 
two parameters if/ and Af are allowed to vary, with the constraint 
that for models of equal mass the same final amount of stars is 
produced. In particular, we constrain the models "H" to produce, 
at the end of the SF period, a mass of newly formed stars which 
is twice as much as the mass of the initial stellar disk (namely 
the fraction = M tt „ ew /M,j is set to be equal to 2). For the sec- 
ond set of models, designated with "L", we set fy = 0.5, namely 
(since the pre-existing stellar disk is more massive for this set 
of models) we allow newly formed stars to be only 50% of the 
pre-existing disk at the end of the SF period. The resulting inte- 
grated luminosity of our models is not far from being constant (a 
constant luminosity has been assumed by MF99). Our reference 
models have Af = 500 Myr and a SFR if/ such that, after Af, the 
required fraction fy is obtained. It is obvious that, keeping the 
starburst duration Af constant, a large starburst luminosity (i.e. 
a large if/) produces a stronger outflow, thus in turn larger ejec- 
tion efficiencies of gas and metals (see MF99). It is however not 
obvious if short, intense starbursts are more effective than long- 
lasting, milder SF episodes in getting rid of gas and processed 
matter. The answer to this question will be sought by consider- 
ing the dependence of the ejection efficiencies in models with 
varying if/ and Af but with constant if/ X Af. In particular, a few 
models have been calculated in which Af is a factor of 10 smaller 
than the reference value (i.e. Af = 50 Myr) and, consequently, 
the SFR is ten times more intense. Our reference value for At of 
500 Myr is justified by observational studies indicating that the 
starbursts in DGs are not as short as previously thought, lasting 
on average a few hundred Myrs (see e.g. McQuinn et al. 2010a 
I2010bl) . 

The initial gas metallicity of our reference sets of model 
galaxies is zero. This is because we do not have a self-consistent 
evolution of the early phases of our model galaxies and we 
can only suppose that the chemical enrichment of gas and stars 
in this early phase will be limited (particularly for the models 
"H"). Moreover, we do not expect it to be a fundamental pa- 
rameter in determining the fate of metals, since the only way 
a different metallicity can affect the development of a galactic 
wind is through the cooling function and the difference between 
the zero-metal cooling curve and that of a modest metallicity 
(Z=10~ 2 Z say) is negligible. We have nevertheless consid- 
ered also models in which the initial metallicity of the stars (and 
gas!) is observationally determined (in particular we make use 
of the M-Z relations observationally determined by Tremonti et 
al. 120041 for the local Universe). 

As done by Fragile et al. (2004) and Rodriguez-Gonzalez et 
al. (1201 II) . another important parameter in our study is the radius 
over which the feedback from dying stars is redistributed. It is to 
expect that, if the "feedback radius" is very small, the density 
of the (metal-rich) ejected material will be so high that radia- 
tive losses can be very significant and can substantially increase 
the probability that the ejecta remain locked inside the galaxy 
(see e.g. Tenorio-Tagle et al. 2007). It is unclear what happens 
if the feedback radius is larger although, as already mentioned, 
the study of Fragile et al. (|2004) suggests a reduction of ejec- 
tion efficiencies (of both gas and starburst matter). Our reference 
feedback radius is Rf = 200 pc but models with Rf = 50 pc and 
Rf = 1000 pc have been also computed. 

We have also considered the possibility (as we did in RH07) 
that the initial distribution of gas is characterized by inhomo- 
geneities. However, for the sake of simplicity, we consider only 
a random perturbation of the initial galactic density distribution. 




0.0 0.5 1.0 1.5 2.0 

R (Kpc) 



Fig. 2. Initial gas distribution for three models (L8R - upper 
panel; L8M - central panel; L8F - lower panel) characterized 
by different degrees of flattening (see Table|2)- The density scale 
(in g cm -3 ) is on the right-hand strip. 

Table 1. Reference parameter values (common to all 18 basic 
models) 



7 (a) 


R F ' b > 


AfW z^W) 





200 


500 8 



Notes. (n> : Initial gas metallicity; <b> : Feedback radius, in parsec; (c) : 
Duration of the star formation period, in Myr; {d) : Redshift of galaxy 
formation 

Namely, once an equilibrium configuration has been obtained, 
we modify the density p(z, j) in the grid (i, j) in this way: 

p(i,j)=p(i,j)-[l + sR [ - 1 , 1] ], (3) 

where s is a small number (corresponding to the largest possible 
amplitude of the perturbation) and is a randomly gener- 

ated number in the range [—1,1]. A more thorough (and self- 
consistent) treatment of inhomogeneities and a detailed study 
of the effect of boundary conditions (gravitational perturbations, 
external pressure, infall of clouds) is deferred to a future paper. 

The central resolution of the simulations is 4 pc for the mod- 
els with 10 8 and 10 9 M of initial baryonic mass and 2 pc for the 
models with 10 7 M of baryons. The resolution decreases out- 
wards with a ratio between adjacent cells of 1.02. The reference 
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values for some key parameters (common to all 18 basic models) 
are summarized in TableQ]whereas the specific parameter values 
distinguishing those models are recalled in Table [2] In order to 
appreciate the distinction between models with different geome- 
tries, Fig.|2]shows the initial gas distribution of the models L8R 
(upper panel), L8M (central panel) and L8F (lower panel). 



4. Results 

4.1. The reference models 

A snapshot of the evolution of all 18 reference models (see Table 
13 is shown in Fig.[3](after an evolutionary time of 100 Myr) and 
in Fig. |4](after an evolutionary time of 200 Myr). The main fea- 
tures of our model results can be noticed from these two figures: 

- As expected, models characterized by a flat initial gas dis- 
tribution (bottom rows) more easily develop a large-scale 
outflow. The roundish models (top rows) show very weak 
signs of outflows even at 200 Myr, irrespective of the bary- 
onic mass. 

- For models characterized by the same geometry, the initial 
baryonic mass clearly plays a very important role, with much 
more prominent outflows in low-mass galaxies (left-most 
columns for each group of panels). This confirms the overall 
trend of gas ejection efficiency vs. mass shown by MF99 for 
instance. 

- If a large-scale outflow is formed, freshly produced metals 
can be easily lost from the galaxy (it is to remind that at an 
evolutionary time of 100 or 200 Myr the SF is still ongoing). 
This is better seen in Fig.[5]where the distribution of oxygen 
(taken as a proxy of metals) is shown for the models with 
small initial gas fractions (models "L"), after 100 Myr. It is 
clear from this figure that in many model galaxies a large 
fraction of oxygen is leaving the parent galaxy. This will be 
quantified in detail below. 

- Outflows in models with high initial gas fraction (models 
"H"; the ones on the left of Figs. [3]and[4| are clearly less 
prominent compared to the corresponding outflows in low- 
gas models. This is due to the fact that the expansion of 
superbubbles is hampered by larger gas pressures in these 
models and to the fact that the SFRs are lower. 

A more quantitative analysis of the results of the models is 
obtained by inspecting the columns 9-12 of Table [2] Here the 
ejected fractions e of gas and oxygen are indicated for each 
model after 200 and 500 Myr. We simply estimate the retained 
gas fractions as the ratios between the mass of pristine gas con- 
tained in 0.5 r V i r at the time 200 (500) Myr and the same mass 
at the beginning of the simulation. To calculate the retained oxy- 
gen fraction we divide instead the mass of oxygen within 0.5 
r vlr by the total amount of oxygen expelled by dying stars until 
200 (500) Myr. The ejection efficiencies (of gas and oxygen) can 
thus be simply estimated as 1 - rf, where rf indicates the cor- 
responding retained fraction. Although our estimate is approxi- 
mate and can be affected for instance by gas (or oxygen) tem- 
porarily leaving the galaxy and being re-accreted later on, it is 
clear that estimates based on the escape velocities (as in MF99) 
are approximate as well (ejected gas does not evolve ballisti- 
cally), thus we keep our simple definition of ejected fractions, 
bearing in mind the inherent uncertainties, thus focusing on dif- 
ferences in the ejected fractions rather than on absolute values. 
A close inspection of the ejected fractions reveal the following 
properties of the models: 



- In roundish models the ejected fractions of oxygen tend to be 
close or slightly smaller than the ejected fractions of pristine 
gas. This is due to the fact that these models do not largely 
depart from their initial spherical symmetries. Under these 
circumstances, most of the metals remain confined inside the 
(almost spherical) superbubbles (the darker central regions 
in the upper rows of Figs. [3]and|4] see also Fig. |5). A large 
fraction of the pristine gas is swept-up in a relatively dense 
shell surrounding the superbubble (the supershell) and, oc- 
casionally, a fraction of this gas can be located at distances 
from the galactic center greater than 0.5 • r Vi >. 

- On the contrary, non-spherical models (models "M" and "F") 
clearly show a tendency of retaining pristine gas more easily 
than metals. This is in agreement with many previous stud- 
ies already mentioned in the Introduction and in Sec. |2](e.g. 
MF99; D'Ercole & Brighenti [T999l Recchi et al. [200Tb . 

- If we focus on models with the same initial gas distribution, 
we recover a clear trend of decreasing metal ejection effi- 
ciencies with increasing galactic baryonic mass, once again 
in agreement with the results of MF99. For instance, the oxy- 
gen ejection efficiencies at the end of the simulations (i.e. 
eo.soo) for the "M" models are shown in Fig. [6] (squares: 
models with low initial gas fractions; triangles: models with 
initial high gas fractions). Oxygen ejection efficiencies are 
close to 1 for models with 10 7 M Q of baryonic mass and 
very low for models with Mb — 10 9 M . 

- As already noticed, models with a larger initial gas fraction 
are less effective in developing galactic winds. Consequently, 
the ejected fractions for these models are systematically 
lower than the corresponding "L" models (the ones with a 
lower initial gas fraction - see also Fig. [5). 

- The effect of geometry on the development of galactic winds 
and on the fate of metals (and also on the fate of pristine gas) 
is quite evident from this table: models with the same bary- 
onic mass can show very different retained fractions depend- 
ing on their degrees of flattening. Model L8R for instance 
retains almost all the oxygen at 200 Myr (eo,200 = 0.001). 
The model L8F, with the same baryonic mass but with a flat 
initial distribution retains only 18.8% of oxygen, i.e. it expels 
8 1 .2% of the produced oxygen through galactic winds. If we 
look at the final ejection efficiencies, we can notice again 
that the roundish model has the smallest oxygen ejection ef- 
ficiency. The models L8F and L8M have both quite low re- 
tained oxygen fractions (therefore high oxygen ejection ef- 
ficiencies), although eo,500 of model L8F is slightly smaller 
than the ejected fraction of L8M. The final retained gas frac- 
tions for all models depend very little on the degree of flat- 
tening. In contrast, focusing on the "H8" family of models 
(high gas fraction, 10 8 M of initial baryonic mass), we can 
notice again that the final ejection efficiencies of gas depend 
very little on b/a, whereas the final oxygen ejection efficien- 
cies are 0.58, 0.375 and 0.061 for the models H8F, H8M and 
H8R, respectively (see column 12 of Tabled, i.e. the metal 
ejection efficiencies can change by up to one order of magni- 
tude depending on the geometry. Similar trends with the de- 
gree of flattening (little effect on the gas ejection efficiencies, 
large effect on the metal ejection efficiencies) can be noticed 
for the families of models "L9", "H9" and "H7", too. The 
comparison of the final oxygen ejection efficiencies for the 
families of models "L9" and "H9" is shown in Fig. [7] 

The last two columns of Table [2] show the oxygen abun- 
dance measured as 12+log(0/H)abundance ratio, where O/H is 
the abundance ratio in number, of all model galaxies. This abun- 
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Table 2. Parameter values for individual models and some representative results. 
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0.6 


0.5 


2.41 


1.43 


1.36 


0.000 


0.980 


0.999 


0.993 


9.84 


9.84 


L7M 


10 7 


1.0 


0.6 


0.5 


2.41 


1.43 


0.85 


0.999 


0.978 


0.999 


0.993 


9.85 


9.81 


L7R 


10 7 


5.0 


0.6 


0.5 


2.41 


1.43 


0.39 


0.993 


0.945 


0.999 


0.991 


8.74 


9.64 


L8F 


10" 


0.2 


0.6 


0.5 


12.4 


2.47 


2.26 


0.365 


0.812 


0.606 


0.826 


7.23 


7.84 


L8M 


10 s 


1.0 


0.6 


0.5 


12.4 


2.47 


1.60 


0.296 


0.111 


0.598 


0.872 


7.99 


7.73 


L8R 


l(r 


5.0 


0.6 


0.5 


12.4 


2.47 


0.90 


0.192 


0.001 


0.558 


0.599 


7.99 


8.21 


L9F 


io 9 


0.2 


0.6 


0.5 


63.5 


4.26 


11.0 


0.041 


0.006 


0.183 


0.639 


7.87 


7.90 


L9M 


io 9 


1.0 


0.6 


0.5 


63.5 


4.26 


5.50 


0.031 


0.001 


0.116 


0.329 


7.88 


8.14 


L9R 


io 9 


5.0 


0.6 


0.5 


63.5 


4.26 


3.45 


0.059 


0.000 


0.123 


0.163 


7.90 


8.26 


H7F 


IO 7 


0.2 


0.9 


2.0 


3.22 


1.58 


1.50 


0.883 


0.938 


0.586 


0.965 


7.11 


6.71 


H7M 


IO 7 


1.0 


0.9 


2.0 


3.22 


1.58 


0.83 


0.860 


0.824 


0.833 


0.943 


7.48 


7.31 


H7R 


IO 7 


5.0 


0.9 


2.0 


3.22 


1.58 


0.47 


0.412 


0.007 


0.627 


0.703 


7.60 


7.68 


H8F 


10 s 


0.2 


0.9 


2.0 


16.5 


2.72 


2.45 


0.057 


0.001 


0.310 


0.580 


7.41 


7.60 


H8M 


10 s 


1.0 


0.9 


2.0 


16.5 


2.72 


1.50 


0.085 


0.000 


0.212 


0.375 


7.42 


7.69 


H8R 


10 s 


5.0 


0.9 


2.0 


16.5 


2.72 


0.94 


0.068 


0.000 


0.170 


0.061 


7.42 


7.84 


H9F 


10 9 


0.2 


0.9 


2.0 


84.7 


4.69 


6.00 


0.025 


0.000 


0.058 


0.086 


7.40 


7.77 


H9M 


10 9 


1.0 


0.9 


2.0 


84.7 


4.69 


3.80 


0.013 


0.000 


0.048 


0.001 


7.39 


7.81 


H9R 


10 9 


5.0 


0.9 


2.0 


84.7 


4.69 


2.45 


0.023 


0.000 


0.042 


0.000 


7.40 


7.81 



Notes. {a) : Initial baryonic mass in M Q (within 0.5 • r vir )\ {h) : Ratio between the scale lengths a and b in the initial Miyamoto-Nagai stellar 
distribution (see Eq.[2}; w : Gas-to-baryon fraction; W) : Ratio between the mass of newly formed stars at the end of the SF period and the mass of 
the pre-existing disk M d ; (e) : Mass of the DM halo (in 10 s M Q ); { - f) : Virial radius (in Kpc, see eq.Q}; (s) - Central gas density (in 10 -24 g cm -3 ); 
Ejected gas fraction after 200 Myr; (,) : Ejected oxygen fraction after 200 Myr; Ejected gas fraction after 500 Myr; Ejected oxygen fraction 
after 500 Myr; (/) : Abundance of oxygen measured as 12+log(0/H), where O/H is the abundance ratio in number, at 200 Myr; (m) : Abundance of 
oxygen measured as 12+log(0/H), where O/H is the abundance ratio in number, at 500 Myr. 
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Fig. 3. Gas density distribution for the 18 reference models (see Table [2) after 100 Myr of evolution. The 9 panels on the left 
represent the models with high initial gas fraction (models "H"); the 9 on the right are the corresponding models with low initial 
fraction of gas (models "L"). For each sequence of 9 panels, the first column represents models with 10 7 M of initial baryonic 
mass; the middle column shows the gas distribution for models with mass 10 8 M and finally the on the right column the models 
with 10 9 M are displayed. The top rows of models are characterized by a roundish initial distribution (models "R"; with b/a = 5). 
The middle rows show models with b/a = 1 (models "M") and finally the bottom rows are characterized by b/a = 0.2 (flat models 
or "F"). At the top-right corner of each panel the model designation is also indicated. For each set of 9 models, the left-hand strip 
shows the (logarithmic) density scale (in g cm -3 ). 



dance ratio is a mass-weighted average of the abundance ratios 
in each computational cell lying within a sphere of radius 0.5 
• r v ,, r . Again, this is only an approximate measure of the galac- 
tic metallicity and only relative abundance ratios between dif- 
ferent model galaxies are relevant. In particular, the present grid 
of models is not aimed at reproducing the observed metallicities 
of individual DGs and it can be noticed e.g., that large galax- 
ies do not always show larger O/H ratios compared to smaller 



galaxies (at variance with the observed M-Z relation; see Lee et 
al. 120061 1. It is also worth noticing that, unexpectedly, the mod- 
els with the highest metallicities are the ones experiencing the 
strongest galactic winds (models "L7"), at variance with the gen- 
eral idea that galactic winds can keep the metallicities of galax- 
ies low (see the Introduction). This is due to the fact that these 
model galaxies expel almost all pristine gas (see Table 13 and 
therefore the (continuously produced) oxygen is mixed with a 
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Fig. 4. Same as Fig. [3]but for snapshots taken after 200 Myr of galactic evolution. 
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Fig. 5. Distribution of oxygen mass density after 100 Myr of 
evolution for the 9 models characterized by a small initial gas 
fraction (models "L" in Table[2]i. The oxygen density scale (in g 
cm 3 ) is on the left-hand strip 



Fig. 7. Oxygen ejection efficiencies after 500 Myr of evolution 
for high-mass (Mb = 10 9 M ) models as a function of the degree 
of flattening b/a. Symbols as in Fig. [6] 
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Fig. 6. Oxygen ejection efficiencies after 500 Myr of evolution 
for medium (b/a = 1 ) models as a function of the initial baryonic 
mass. Squares refer to models with initially low gas fractions and 
triangles to models with high f g . 



very small fraction of unpolluted gas. The average metallicity 
of the galaxy is thus close to the (very large) metallicity of the 
stellar ejecta. However, this gas is deemed at leaving the par- 
ent galaxy and is anyway too hot to form star; it will thus not 
increase the average metallicity of the stellar populations in the 
galaxy. The models attaining the lowest metallicities are the ones 
(like H7F) developing large-scale outflows but retaining a large 
fraction of pristine gas bound to the galaxy. See also Sect. |5]for a 
more extended discussion on the connection between the results 
of our models and the M-Z relation. 



4.2. Widening the parameter space 



As already mentioned in Sect. 13.21 we have considered also mod- 
els with different SFHs, different feedback radii, different initial 
metallicities and with perturbed initial distributions of gas. We 
summarize in this subsections how these parameters affect our 
findings. 
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4.2.1. Dependence on the star formation history 

We have re-run some models (specifically the ones for the family 
"L8") in which the SFH, instead of being constant for 500 Myr 
as in the standard models, is different from zero only during the 
first 50 Myr, but the SFR is ten times more intense than in the 
standard models. A large-scale outflow quickly develops, even 
for the roundish model, sustained by the large rate of SNII explo- 
sions during the first ~ 80 Myr. At the end of this phase, energy 
is still supplied by Type la SNe and winds from intermediate- 
mass stars, but at a lower rate. This energy is still able to sustain 
the outflow for some time, but after a few times 10 8 yr, the grav- 
itational pull prevails over the pressure gradient caused by the 
hot cavity of gas and some gas begins to recollapse towards the 
center of the galaxy (see Recchi & Hensler I2006I for a precise 
assessment of this phenomenon). 

The result is that the calculated ejection efficiencies (of gas 
and metals) are larger at 100 or 200 Myr than at 500 Myr. In par- 
ticular, the ejection efficiencies at 200 Myr are much larger than 
the ones in the standard models, tabulated in Table [2] Disk-like 
models (models "F" and "M") retain only a few percent of gas 
and metals, and also the roundish model ("R") is able to expel 
more than 50% of the initial gas at distances larger than 0.5 • r V! >. 
The large gas ejection efficiencies for these models are due to the 
fact that, at these high rates of energy release, also the horizontal 
transport of gas can be effective and some gas is pushed to large 
distances from the galactic center even along the disk. The disk- 
like "M" and "F" models show again metal ejection efficiencies 
(slightly) larger than gas ejection efficiencies. Consequently, the 
metallicities of these model galaxies are quite low during the SF 
episode. 

However, as soon as the luminosity considerably reduces, 
this gas is pulled again within 0.5 ■ r V! > (i.e. a fall back is oc- 
curring), whereas the material previously channelled along the 
galactic wind accelerates away of the galaxy due to the steep 
density gradient. Consequently, after 500 Myr of evolution, the 
retained gas fractions for these models tend to be larger (by a fac- 
tor of ~ 1.5-2) than the corresponding fractions of the standard 
models (the ones reported in column 1 1 of Tabled, whereas the 
retained fraction of metals are only slightly larger. 

4.2.2. Dependence on the feedback radius 

We have considered (again for the family of models "L8") a vari- 
ation of the radius over which energy and metals are injected. As 
already anticipated in Sect. 13.21 instead of the reference value 
Rf — 200 pc, we have calculated models with Rf = 50 pc and 
Rf = 1000 pc. The models with Rf = 1000 show very large 
ejection efficiencies for both gas and metals. This is due to the 
unfortunate circumstance that 0.5 • r Vi > for this family of models 
is only slightly larger than 1000 pc, therefore a moderate en- 
ergy injection already suffices to push a large fraction of gas at 
distances larger than 0.5 • r w >. The comparison with the ejec- 
tion efficiencies of standard models makes thus little sense. We 
have therefore calculated the ejection efficiencies of these mod- 
els (and of the corresponding standard models) at r — r,,,>. These 
turn out to be ~ 10% smaller than the ones of corresponding 
reference models. We confirm therefore that, if the energy is re- 
distributed over a large volume, radiative energy losses are more 
effective and the final ejection efficiencies are reduced, in agree- 
ment with the findings of Fragile et al. (120041 1. Also the mod- 
els with Rf = 50 pc show slightly lower ejection efficiencies. 
In this case, the enhanced cooling efficiency is due to the fact 
that the gas within the feedback radius is now characterized by 



large densities and large metallicities (see also Tenorio-Tagle et 
al. 120071 1. Also in this case, the effect of the feedback radius is 
quite limited (ejection efficiencies change by ~ 10-15 %). As 
also shown by Rodriguez-Gonzalez et al. (1201 II see their fig. 7) 
the dependence of the ejection efficiencies on the feedback ra- 
dius is non-monotonic, although the dependence of our results 
on this particular parameter is quite limited. 

4.2.3. Dependence on the initial metallicity 

Since we consider models in which a stellar disk is already 
present at the beginning of the simulations, it seems unreason- 
able to start with a primordial metallicity. We considered there- 
fore models in which the initial metallicity (of stars and gas) is 
regulated by the mass of the pre-existing disk, namely, according 
to the M-Z relation, the larger the disk mass, the larger the initial 
metallicity we must consider. In particular, if we take into ac- 
count the correlation between mass and metallicity obtained by 
Tremonti et al. d2004l l. models of the "L8" family should have 
an initial metallicity of about one tenth of solar. Of course, the 
final metallicities of these models will be much higher than the 
ones attained by models initially without metals. However, as al- 
ready mentioned, the only dynamical effect of a different initial 
metallicity is to increase the radiative losses and thus to reduce 
the amount of thermal energy of the galaxy. It turns out that this 
reduction is quite limited and, therefore, the ejection efficiencies 
do not change substantially (they are only a few per cent lower 
than the ejection efficiencies of corresponding models initially 
without metals). For these models, the metal ejection efficiencies 
are calculated according only to the amount of freshly produced 
metals retained or ejected. 

4.2.4. Models with random perturbation of the initial gas 
distribution 

As described in Sect. 13.21 once an equilibrium initial configura- 
tion has been obtained, we have perturbed it by means of eq.|3] In 
particular, we have perturbed the model L8R by amplitudes s of 
1% (mildly perturbed model) and 5% (largely perturbed model). 
The perturbation increases the turbulence of the model. Mixing 
between hot and cold regions is enhanced and the net effect is an 
increase in the radiative energy losses, hence a reduced amount 
of energy available to drive galactic winds (in other words, the 
galactic winds become more mass loaded). As a consequence 
of that, the retained gas fraction of the mildly perturbed model 
increases (from the final value of 0.442 of the model L8R it in- 
creases to 0.547). However, the final retained fraction of freshly 
produced metals (slightly) decreases for this model (it reduces 
to 0.345 from the reference value of 0.401). This is due to the 
fact that it is much more difficult to keep the model close to 
spherical symmetry if the initial gas distribution is perturbed. 
Regions of lower pressure can be created randomly within the 
computational box and the propagation of metals can thus de- 
viate from isotropicity, being faster (and allowing some venting 
out of metals) along (randomly oriented) directions with steep 
pressure gradients. This effect is similar to the one analyzed in 
great detail by RH07 (see also Sect. [2] for a short summary of 
the main results of this paper). In the largely perturbed model, 
the final metal ejection efficiency is similar to the one obtained 
in the mildly perturbed model. However, the venting out (along 
random directions) involves now also a non-negligible fraction 
of pristine gas and the final gas retained fraction is 0.443 (very 
similar to the gas retained fraction of the reference model L8R). 
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Fig. 8. Distribution of gas (left panel) and oxygen (right panel) 
density for a model with large perturbation of the initial den- 
sity distribution (see text), after 200 Myr of evolution. Brighter 
colors correspond to larger densities. 



The piercing of the supershell along different directions and not 
just along the direction perpendicular to the disk of the galaxy 
can be appreciated by inspecting Fig. [8] It shows the gas and 
oxygen density distribution after 200 Myr of evolution of the 
largely perturbed model. Clearly, the perturbation of the initial 
gaseous distribution affects the ejection efficiencies of gas and 
metals in a non-linear (and some times non-predictable) way and 
deserves further studies. 

It is important to stress that, due to the assumed symmetry 
and dimensionality of the code, inhomogeneities represent rings 
of denser gas. It is not clear whether a realistic 3-D distribu- 
tion of inhomogeneities will produce the same results. Our group 
is currently running (with the AMR code FLASH) 3-D chemo- 
dynamical simulations of DGs with inhomogeneities (Mitchell 
et al., in preparation). The results of these simulations can shed 
light on the reliability of our results. 

5. Discussion and conclusions 

In this paper we have studied the development of galactic out- 
flows in DG models. Our main focus is to study the dependence 
of the ejected mass fractions (of gas and freshly produced met- 
als) on the degree of flattening of the galaxy, a problem that has 
received little attention in the past. A very simple theoretical 
analysis leads to the conclusion that, in a very flat galaxy, the 
ejection efficiencies of freshly produced metals must be quite 
high. In fact, once the over-pressurized gas formed by SN ex- 
plosions and stellar wind breaks out and forms a bipolar galac- 
tic wind, metals can be easily ejected out of the galaxy through 
this funnel. On the other hand, in a spherical (or almost spher- 
ical) galaxy, bipolar outflows cannot be formed because there 
is no preferential direction, along which the pressure gradient is 
steeper. Therefore, either the galaxy expels gas isotropic ally (but 
blow-away is very difficult in DGs, see e.g. Hensler et al. 2004) 
or the superbubble of hot (and metal-enriched) gas remains con- 
fined inside the galaxy (thus the metal ejection efficiencies are 
very low). We have considered three different degrees of flatten- 
ing for our model galaxies: a roundish model (one for which the 
ratio between the length scales a and b describing the Miyamoto- 
Nagai potential is b/a - 5), a thick-disk model (with b/a = 1) 
and a thin disk model (with b/a = 0.2). Our study confirms the 
trend described above: ejection efficiencies in roundish galaxies 
are systematically lower than the ones in disk-like galaxies (see 
Table [2] or Fig.|7]i. On the other hand, for most of the analyzed 
models, transport of gas along the disk is quite limited, therefore 
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Fig. 9. Oxygen abundances (in number) after 500 Myr of evolu- 
tion for high gas fraction models (models "H", with f g = 0.9), as 
a function of initial baryonic mass. Squares refer to flat models 
(models with b/a = 0.2); triangles show the results for inter- 
mediate models (b/a = 1) and rhombs correspond to roundish 
models (b/a = 5). 



(in agreement with many previous studies, e.g. MF99; D'Ercole 
& Brighenti 1999) we can conclude that ejection efficiencies of 
gas are (at least for disk-like models) lower than ejection effi- 
ciencies of freshly produced metals. 

We have also studied in detail the effect of galactic mass 
on the fate of gas and freshly produced metals, a problem 
that, on the contrary, has received much more attention in the 
past (MF99; Rodriguez-Gonzalez et al. 1201 11 1. We confirm that 
smaller DGs, with shallower potential wells, favor the develop- 
ment of large-scale outflows, hence the ejection efficiencies (of 
gas and metals) increase with decreasing galactic masses (see 
Table [2] or Fig. |6). The increase of metal ejection efficiencies 
through galactic winds has been often invoked to explain the ob- 
served M-Z relation in galaxies (see e.g. Tremonti et al. 2004). 
The dependence of the obtained metallicity (measured by means 
of 12+log (O/H)) of our model galaxies as a function of their 
initial masses is shown in Fig.[9]for the "H" set of models. For 
all three degrees of flattening, a trend of increasing metallicity as 
a function of mass is clearly visible. However, as outlined also in 
Sec. 14. II for the "L" set of models this trend is much less visible. 

It is important to notice that the abundances tabulated in 
Tableland plotted in Fig.|9]are mass-averaged means within the 
whole galactic region. However, two gas phases coexist within 
this region: a hot gas phase (the galactic wind and the starburst 
region) and a colder ISM. We expect the hot phase to have much 
higher metallicities than the cold phase. The metallicity of the 
galactic wind is expected to play no role in the process of chem- 
ical enrichment of the galaxy (this gas is destined to leave the 
parent galaxy) but is important in the process of chemical enrich- 
ment of the infra-cluster medium. To have an idea on the oxygen 
abundances of these two phases, we distinguish grid points with 
temperatures above and below 2 • 10 4 K. Gas above this tem- 
perature threshold is supposed to be too hot to be detected by 
optical spectroscopy. The resulting values of 12+log(0/H) after 
500 Myr of evolution for the "H" family of models is shown 
in Fig. [10] (red lines: cold phase; blue lines: hot phase). As ex- 
pected, the abundances in the cold phase are much lower than the 
abundances in the hot phase. The cold phase abundances show 
the same trend of the total abundances shown in Fig. [9] This is 
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Fig. 10. Oxygen abundances (in number) after 500 Myr of evo- 
lution for "H", as a function of initial baryonic mass. Symbols 
are as in Fig. [9] The abundances in the cold (T < 210 4 K) phase 
are shown in the lower set of points (connected by red lines). 
The abundances in the hot (T > 2-10 4 K) phase are shown in the 
upper set of points (connected by blue lines). 



mainly due to the fact that the calculated mean abundances are 
mass-weighted and most of the gas mass is colder than 2 ■ 10 4 
K. For some models (in particular the models with log M=7), the 
cold phase abundances are significantly lower (~ 0.2 dex) than 
the total abundances. This is expected because, in some models, 
the metal-rich hot gas (wind and starburst region) represents a 
significant fraction of the total gas. On the other hand, models 
experiencing no (or very limited) galactic winds contain much 
less hot gas, therefore the average total abundance is less af- 
fected by the metallicity of the hot phase and is much closer 
to the abundance of the cold gas (the differences are a few hun- 
dredths of dex). 

Also the abundances of the hot phase tend to grow with mass 
(see the blue lines in Fig. ITOb . High mass galaxies show higher 
hot phase abundances because, for these models, much of the 
hot gas is in the starburst region and thus its metallicity is close 
to the (very high) metallicity of the stellar ejecta. On the other 
hand, galaxy models experiencing strong galactic winds show 
smaller hot-phase metallicities because the winds entrain more 
cold (and metal-poor) gas. This process is usually indicated as 
mass-loading and is confirmed by many studies (both theoreti- 
cal and observational) of galactic winds (see e.g. Strickland & 
Steve nsEOOOl Tescari et al. 120091 Hopkins et al. l20T2l Newman 
et al. 120121 among many others). Indeed, a clearer correlation 
can be shown between the ejected oxygen fraction after 500 Myr 
(12th column in Table |2]i and the hot-phase oxygen abundance. 
This correlation is shown in Fig.Q~T] Models experiencing galac- 
tic winds (thus with high oxygen ejection efficiencies) show also 
mass-loading, thus dilution of the oxygen abundances in the hot 
phase. On the other hand, in models with weak (or without) 
galactic winds most of the hot gas is concentrated in the star- 
burst region, where the abundances closely resemble the (very 
high) abundances of the stellar ejecta. 

One should warn the reader that the processes of metal dis- 
persal and mixing are very difficult to model. No hydrodynam- 
ical code is able to correctly resolve contact discontinuities that 
form between hot cavities and cold shells. Numerical diffusion 
tends to smear this discontinuity and create regions of interme- 
diate temperature and densities, where metals can cool in a rela- 
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Fig. 11. Oxygen abundances in the hot phase plotted versus the 
oxygen ejection efficiencies (12th column in Table [2). Model 
designations are labelled close to each point. 



tively short time-scale and then mix with the surrounding shell. 
In spite of these uncertainties, we believe that 500 Myr is a suf- 
ficient time span to cool down most of the produced metals (if 
they are not channelled in a galactic wind). This is confirmed 
by some of our previous studies (see in particular Recchi et al. 
120011 their sect. 3.2) but also by other studies addressing the 
specific issue of metal mixing (e.g. de Avillez & MacLow 2002; 
Pan et al. 2012; Yang & Krumholz 2012), which confirm that a 
few hundred Myrs is a sufficient time to disperse and mix metal- 
licity inhomogeneities. It is also important to stress that, in our 
results, cooling and mixing of stellar ejecta is mainly due to ther- 
mal conduction and to the formation of eddies and vortices. The 
presence of inhomogeneities (see Sect. 14.2.4b can facilitate both 
processes, thus enhance the mass loading of the galactic winds. 

One should also notice that the chosen temperature thresh- 
old separating the hot and the cold phase (2 • 10 4 K) is some- 
what arbitrary and it has been taken in analogy with what done 
in some of our previous publications (see in particular Recchi et 
al. 2006). A different and maybe better choice could be 4 • 10 4 K 
(the minimum of the cooling curve for low values of the metal- 
licity - see Bohringer & Hensler[1989 their fig. 1), or even 10 5 
K. We have checked how much are our results affected by the 
choice of the temperature threshold. It turns out that the abun- 
dances in the cold phase are almost unaffected by this choice 
(for T = 4 • 10 4 K the abundances change negligibly whereas for 
T = ■ 10 5 K the differences are 0.02 dex at most). This is due to 
the fact that a small fraction (in mass) of gas has temperatures 
in the range [2 • 10 4 , 4 • 10 4 ] K or [2 ■ 10 4 , 10 5 ] K. On the other 
hand, the variations in the hot phase abundances are more sig- 
nificant. They increase by up to 0.15-0.2 dex if the threshold is 
10 5 K. However, the overall correlations shown in Figs. [TOland 
fTTIremain. 

We must outline once again that the main aim of this paper 
was to show the effect that the gas distribution can have on the 
development of galactic winds and, hence, on the fate of gas and 
freshly produced metals. We have not attempted to reproduce 
the observed M-Z relation (we will focus on this task in a pa- 
per in preparation). For instance, all the models, irrespective of 
their masses, have the same SFRs per unit available gas mass 
(hence the same SF efficiencies). Models with initially 10 8 M Q 
of baryonic mass have SFRs ten times larger than models with 
Mb = 10 7 M Q . The metal production rates are also ten times 
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larger. However, these larger fractions of metals are mixed with 
ten times more pristine gas. Hence, if galactic winds do not play 
a significant role, the final metallicities of the models should be 
approximately independent on the galactic mass (see roundish 
models in Fig . |9j> ■ However, it is known from chemical evolution 
studies (e.g. Matteucci fl994t Pipino & Matteucci|2004]) that SF 
efficiencies must increase with galactic masses in order to re- 
produce the observed chemical abundances and abundance ra- 
tios of galaxies. This is also consistent with the idea of having 
SF efficiencies increasing with pressure of the ambient diffuse 
gas, which increases with galactic mass (Elmegreen & Efremov 
[19971 Harfst et al. l2006l see also Leroy et al. 120081 

However, we have noticed in Sect. 14.2.11 that metal ejec- 
tion efficiencies of galaxies with higher SFRs (but producing at 
the end of the simulation the same amount of stars) are larger. 
Consequently, these galaxies tend to have lower metallicities 
during the SF episode. This qualitative trend is in agreement with 
the so-called fundamental metallicity relation (FMR: Mannucci 
et al. 1201 Oi l, according to which the oxygen abundance in star 
forming galaxies correlates with the quantity fi a = log(M„) - 
a \og(S FR), where a is a free parameter chosen to minimize the 
scatter in the FMR. Although the precise shape of the FMR is de- 
bated (for instance the recent paper of Andrews & Martini BO 121 
find a = 0.66, more than twice the value of Mannucci et al.), un- 
doubtedly the SFR plays a very important role in determining the 
gas-phase metallicity of a galaxy. Galaxies with the same stellar 
mass but with higher SFRs are characterized by lower metallic- 
ities (see e.g. fig. 1 1 of Andrews & Martini 2012 1. According to 
Andrews & Martini (2012), this is probably due to the fact that 
galaxies with high SFRs are presently experiencing a merging. 
Major mergers drive in considerable amounts of low metallicity 
gas from large radii, which dilutes the metallicity of the galaxy 
and triggers a vigorous SF burst. Since most of the galaxies with 
large SFR have already large stellar masses, it seems unlikely 
that their gas components have still a low metallicity. It seems 
to us more likely that large SFRs drive large outflows, with a 
consequent significant metal loss, as our models show. 

What we want to point out here is that, for models with the 
same mass and SFR, geometry plays a very significant role in de- 
termining the fate of freshly produced metals and, consequently, 
the final metallicity. An inspection of Table [2] shows that, for 
some models, ejection efficiencies of metals can change by up 
to an order of magnitude depending on the degree of flatten- 
ing, being instead the gas ejection efficiencies quite independent 
on b/a. Consequently, the final metallicities of models with the 
same mass can vary up to 1 dex depending on the geometry (see 
Figs. l9l and fTOb. The large spread observed in the metallicity of 
DGs with the same masses (see e.g. Lee et al. 120061 their fig. 
8) could be due to the effect of gas distribution. Although the 
parameter a is chosen to minimize the scatter in the FMR, still 
galaxies with the same stellar mass and SFR show a spread in 
the metal content (see e.g. figs. 1 1 and 12 of Andrews & Martini 
120121 1. The gas distribution might be responsible for this spread. 

A last comment related to the M-Z relation (or to the FMR) 
concerns the correlation between mass and degree of flattening 
of DGs. Recent studies suggest that smaller DGs tend to have 
larger axial ratios (i.e. they tend to be rounder) than DGs with 
larger masses. This might complicate the interpretation of the 
M-Z and FMR relations according to our models. However, this 
effect seems to be quite limited. For instance, Sanchez-Janssen et 
al. (1201 Oi l show (their fig. 1) that the average axial ratios change 
by at most ~ 20% between DGs with stellar masses of 10 7 and 
10 9 M , and that the spread in axial ratios is extremely large at 
all stellar masses. Moreover, Lisker et al. (2009 i show that (at 



least for DGs in clusters) the axial ratio is more related to the 
DG velocity than to the mass. 

Our main conclusions can be summarized as follows: 

- The gas distribution in a galaxy plays a very important role 
in determining the fate of freshly produced metals. Model 
galaxies with the same masses but with different degrees of 
flattening can have metal ejection efficiencies differing for 
one order of magnitude. In particular, flat galaxies easily 
develop bipolar outflows, through which a large fraction of 
metals can be lost. 

- On the other hand, the fate of pristine gas is less depen- 
dent on geometry than the fate of metals. Since some models 
(with the same initial mass) can have similar ejection effi- 
ciencies of pristine gas but very different ejection efficiencies 
of metals, the final attained metallicities can vary by up to 
one dex (being the roundish models characterized by larger 
final metallicities). 

- Models characterized by the same degree of flattening show 
a clear dependence of the metal ejection efficiencies on the 
galactic mass. Smaller galaxies (with shallower potential 
wells) more easily develop large-scale outflows, therefore 
the fraction of lost metals tends to be higher. 

- Ejection efficiencies (of gas and freshly produced metals) 
significantly depend also on the star formation history of the 
galaxy and on the presence (or absence) of density perturba- 
tions. 

- The ejection efficiencies show instead a moderate depen- 
dence on the initial metallicity and on the size of the region in 
which the energetic and chemical feedback from Supernovae 
and stellar winds is redistributed. 

Acknowledgements. We thank the anonymous referee for very helpful com- 
ments and remarks. 
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